Non-sequential double ionization of molecules 
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Double ionization of diatomic molecules by short linearly polarized laser pulses is analyzed. We 
consider the final stage of the ionization process, that is the decay of a highly excited two electron 
molecule, which is formed after re-scattering. The saddles of the effective adiabatic potential energy 
close to which simultaneous escape of electrons takes place are identified. Numerical simulations 
of the ionization of molecules show that the process can be dominated by either sequential or non- 
sequential events. In order to increase the ratio of non-sequential to sequential ionizations very short 
laser pulses should be applied. 
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I. INTRODUCTION 

High intensity ultrashort-pulse lasers allow experi- 
mental studies of multi-electron effects, such as a non- 
sequential double ionization, high order harmonic genera- 
tion or above threshold ionization [JQ- While the single 
ionization of atoms or molecules, as well as the high order 
harmonic generation can be described within a single ac- 
tive electron model, such an approximation in the case of 
double ionization and laser intensity below the saturation 
value gives ionization rates that are much smaller than 
experimentally observed, indicatingthat interactions be- 
tween electrons are important P,0iH,0- In this paper we 
consider the double ionization of molecules within a clas- 
sical model for electrons in a combined Coulomb and ex- 
ternal field applying the approach developed in 00,00] 
for the multiple ionization of atoms. 

Multi-photon double ionization of atoms in strong laser 
fields can be regarded as a three step process 0, Q • In 
the first step one electron tunnels out through the Stark 
saddle and then is returned back to the nucleus 0, ^| . 
Thereafter, in the second step a highly excited state of 
an atom or a molecule is formed at the expense of the 
energy brought back by the returning electron (up to 
i.nUp, where U p is a ponderomotive energy). Finally, in 
the third step such a highly excited compound state de- 
cays in several ways through a single, double or multiple 
ionization. The starting point in our classical analysis of 
the double ionization of molecules is the excited complex, 
i.e. we assume that we have an initial state of two highly 
excited electrons close to the molecular core. We shall fo- 
cus on double ionization events that can then appear as a 
possible channels of decay. Among those events there are: 
a sequential process (when the highly excited compound 
state decays emitting electrons one by one) and a non- 
sequential event (two electron leave the excited molecule 
simultaneously) . 

There have been experimental research aimed on dou- 
ble ionization of diatomic molecules 0, 0, 0, Q 
which showed that there are differences between molecu- 
lar species. For example, nitrogen molecule, N2, clearly 



exhibits a "knee-structure" in a double-to-single ioniza- 
tion yields ratio, like atoms do, while for oxygen molecule, 
O2, such a structure is rather not visible |12.[l3| . More- 
over, in the case of N2 it seems, that electrons escape 
with similar momenta along field polarization axis more 
often than in the case of O2 01 • We consider this prob- 
lem from the point of view of the classical analysis, which 
starts after the formation of the highly excited compound 
state, and discuss possible reasons for such different ex- 
perimental observations. 

Paper is structured as follows in Sec. [H] we intro- 
duce a theoretical model on which we based our anal- 
ysis. Sec. IIIII consists of the identification and descrip- 
tion of the saddles in the potential. Then in Sec. IIVI we 
present results of the numerical simulations. And finally, 
in Sec. we conclude. 



II. MODEL 

As we mentioned above, in our starting situation (as 
a result of the former re-scattering) there are two highly 
excited electrons close to the molecular core in the pres- 
ence of linearly polarized laser field. The motion of the 
molecular core is frozen due to the fact that for short 
laser pulses molecules have not enough time to change 
their orientation |14l Il5|. Thus, the Hamiltonian reads 
(in atomic units): 



where the potential, 



V = V 1 + V 2 + V 12 + V F , 



(1) 



(2) 



consists of the potential energies associated with interac- 
tions of the electrons with the nuclei (the entire structure 
of the molecular core is approximated by two positively 
charge nuclei), 
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for (i = 1,2), where R; and indicate position of the nu- 
clei and the electrons, respectively, the repulsion between 
electrons, 
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and the term describing interaction with the field (polar- 
ized along the z axis), 



V F =F{t)(z 1 + z 2 ). 



(5) 



The electric field strength F{t) has an oscillatory com- 
ponent times the envelope from the pulse, namely: 



F(t) = Ff(t) cos(wt + . 



(6) 



with F, lo and <p the peak amplitude, frequency, and 
initial phase of the filed, respectively, and with 



f(t) =sin 2 (7rf/T d ), 

the pulse envelope of duration 

2vr 

T d = n — , 



(J) 
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where n is number of cycles in the pulse. 

For a sake of convenience, we shall place the origin of 
our coordinate system in the center of mass of the nu- 
clei, then two parameters appear, namely: d, distance 
between nuclei, and angle between molecular axis and 
z axis (polarization axis). Without loss of generality we 
assumed that the molecule lies in xz plane, hence poten- 
tial energies for each electron is: 
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As one can see, in our model the only difference be- 
tween different diatomic species lies in the distance, d, 
between the nuclei. 



III. LOCAL ANALYSIS 



atoms, those two saddles lie symmetrically with respect 
to the polarization axis. Then the motion of the elec- 
trons can be confined in some symmetry subspace 0, 0| • 
For diatomic molecules (such as N 2 or 2 ) that does not 
generally occur since they possess their own symmetry 
axis which can be oriented at any angle to the polariza- 
tion axis destroying the global symmetry. Nevertheless, 
saddles are formed by the external field and following ar- 
guments presented in [HIS 00 we assume that electrons 
to leave the molecule in a correlated manner have to pass 
simultaneously close to them. 

Before the double ionization escape both electrons pass 
close to the nuclei where they interact strongly with each 
other and with the nuclei. For that reason we assume 
that all memory of the previous motion is lost. Then it is 
correct to assume that the compound state which decays 
to a doubly charged molecule may be classically simu- 
lated by a statistical distribution for two electrons close 
to the nuclei. Furthermore, the classical motion of the 
electrons is fast compared to the field oscillations and an 
adiabatic approximation, keeping the field fixed, becomes 
useful in the analysis of the ionization channels. Using 
this adiabatic assumption we will identify and describe 
saddle points, through which the electrons can escape. 

For a molecule oriented along the field axis the problem 
possesses axial symmetry. Then switching to the cylin- 
drical coordinates (pi, tpi, Zi for i— th electron) one can 
easily define a symmetry subspace of electron motion. 
Restricting the electrons to a plane (i.e. </?i — if2 = tt) 
their coordinates in the C 2v symmetry subspace are 
pi = P2 = R, Z\ — Z2 = Z and the potential energy 
reduces to 
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With a molecule oriented perpendicularly with respect 
to the field there are two C 2v symmetry subspaces. One 
subspace is defined in the xz plane the other in the yz 
plane. In the former case, the electron coordinates in 
the subspace are (x% — X, yi — 0, z\ — Z) and (x 2 = 
—X, 2/2 = 0, Z2 = Z) and the potential energy reads 



What we consider in this paper is the evolution of elec- 
trons in combined Coulomb and external fields after a re- 
scattering event when the highly excited compound state 
is formed. For one electron in the Coulomb potential, if 
the external field is non-zero the Stark saddle is opened. 
Then the electron can get away through this saddle and 
ionize. In our case we have two electrons and if there 
were no interaction between the electrons they could es- 
cape simultaneously through the same saddle on top of 
each other. But once there is repulsion between electrons, 
the Stark saddle splits into two saddles that lie on the op- 
posite sides of the field polarization axis. In the case of 



V 



I) 2 



z 2 



(11) 



(X + 1) 2 



z 2 



2\X\ 



2F(t)Z. 



And for saddles which are in the yz plane, namely coor- 
dinates are {x\ — 0, y\ = Y, z\ = Z) and {x% = 0, y 2 = 
—Y, Z2 = Z) and the potential energy is 



V= - 
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The potential energies Eqs. I|1U[1. and i|12|l are shown 
in Fig. n f° r a se t of parameters corresponding to the 
nitrogen molecule (internuclear distance is d = 2.067 
a.u.) and for the field F(t) = 0.07 a.u. (an intensity 
of 1.7 x 10 14 W/cm 2 ). The saddles are clearly visible. 

In a non-sequential double ionization the electrons 
have to pass sufficiently close to the saddles (exempli- 
fied in Fig. ^ for molecular orientation being 8 = and 
8 = 7r/2). For a general orientation of a molecule there 
are two possible channels for non-sequential ionization 
— the electrons escape by passing saddles situated in a 
plane defined by the field and molecular axes or in the 
perpendicular plane. These two channels are equivalent 
when a molecule is parallel to the field and the axial sym- 
metry is restored. 

For a fixed external field local analysis of the saddles 
reveals a few stable and unstable directions. The latter 
are responsible for either simultaneous double electron 
escape or single ionization (in such a process one elec- 
tron leaves the molecule while the other returns to the 
core). Each of unstable directions can be characterized 
by a Lyapunov exponent. Knowledge of the Lyapunov 
exponents allows one, similarly like in the problem of 
double ionization without the external field considered 
many years ago by Wannier 0,0,0,0], to derive the 
dependence of the cross section on energy close to the 
threshold, namely, 



a(E) cx {E-VsT, 
where Vs is a saddle energy and the exponent 



(13) 



(14) 



where A r is the Lyapunov exponent of the unstable di- 
rection corresponding to the non-sequential double ion- 
ization and Xi are all other Lyapunov exponents of a 
saddle 0,0. 

We examine properties of the saddles as the internu- 
clear distance, d, increases — the results are shown in 
Fig. Starting with d = (in that case 8 is meaning- 
less since a molecule reduces to an atom) and increasing d 
the energy Vs of the saddle corresponding to the molecule 
orientation 8 = is always the lowest and it decreases. 
For the other extremal orientation, i.e. 8 = 7r/2, the en- 
ergies of the saddles (there are two saddles because of 
the axial symmetry breaking) increase and their values 
are the highest. Analyzing the dependence of a on the 
internuclear distance we see that increasing d the cross 
section exponent of one of the saddles (corresponding to 
8 = 7r/2) goes up while the other goes down. The expo- 
nent for 8 = slightly increases. 

Now let us examine how the parameters of the saddles 
change with the orientation 8 in the case of N 2 (d = 2.067 
a.u.) and O2 (d = 2.28 a.u.) molecules, see Fig. |3] The 
energy Vs of all saddles increases and the exponent a of 
one member of the saddle pairs increases while the other 
decreases. Fig. [3] shows that the parameters for N2 and 




FIG. 1: Section through equipotent surfaces of the adiabatic 
potential, V = V1 + V2 + V12 + Vf, for fixed time t and for two 
symmetric orientation of the molecule, namely parallel (9 = 0, 
top panel corresponding to Eq. 1101 1 and perpendicular to 
the field polarization axis, Z, (9 = it/2, middle and bottom 
panels corresponding to Eq. li lt and Eq. 1121 . respectively). 
F{t) = 0.07 a.u., d = 2.067 a.u. 



O2 molecules are quite similar — the largest differences 
between the species are of the order of few percent. Tak- 
ing into account that the experimental results 0, 0] 
are the statistical mixtures of different molecule orienta- 
tions we may conclude that from the point of view of our 
local analysis of the non-sequential decay channels there 
should be no differences between N2 and 02- 
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FIG. 2: Energy of the saddle (top panel) and cross section 
exponent (bottom panel) as a function of an increasing inter- 
nuclear distance, d, for two different orientation angles, 8, i.e. 
9 = and 6 — n/2 (the dashed line correspods to the saddle 
in xz plane and the dash-dotted line correspods to the saddle 
in yz plane, respectively). F(t) = 0.07 a.u. 



FIG. 3: Energy of the saddle (top panel) and cross section 
exponent (bottom panel) as a function of an orientation angle, 
9, for N 2 (d = 2.067 a.u.) and 2 (d = 2.28 a.u.) molecules. 
F(t) = 0.07 a.u.. (For both molecules the top line corresponds 
to the saddle in the yz plane, whereas the bottom line to the 
saddle in the xz plane) 



IV. NUMERICAL SIMULATIONS 

So far we have discussed the local analysis of the poten- 
tial within the adiabatic assumption. Now we are ready 
to go one step further and perform numerical simulations. 
In the Coulomb potentials associated with interaction be- 
tween the electrons and nuclei, Eq. (JOJ, we introduce a 
smoothing factor e [2lll2^ | to avoid divergence in numeri- 
cal integration of equation of motions. Then the potential 
terms read 




We choose e = 0.01 which introduces negligible changes 
of the Vs and a parameters presented in Figs. [3 and 

Assuming an initial re-scattering event took place, two 
excited electrons pass close to the nuclei where they inter- 
act strongly with each other and with the nuclei. There- 
fore, it is reasonable to assume, as mentioned in Sec. IIIII 



that all the memory of the earlier motion is lost. Then 
the initial state of the final stage for double ionization is 
a statistical distribution for two electrons close to the nu- 
clei. Hence, we choose initial values of positions and mo- 
menta with respect to the microcanonical distribution for 
a given initial energy E (it should be in the range between 
—I and —I + 3.17£/p, where / is the ionization energy). 
The positions are chosen microcanonically but with ad- 
ditional conditions, namely z-, = and \J :rj + t/f < 85 
a.u. p • We start all simulations at the peak of the laser 
pulse with the phase of the field chosen randomly. In all 
simulations an ensamble of 10 5 trajectories is used. 

Distributions of the final electron momenta parallel to 
the field axis for different initial energy E and for different 
pulse durations are presented in Fig. 0] . All data in the 
figure correspond to the N2 molecule oriented along the 
field axis. For very short laser pulses (n = 2 cycles) sig- 
natures of simultaneous electron escape are clearly visible 
— the distributions are localized along the diagonals in- 
dicating that the electrons escape predominately by pass- 
ing close to the saddles analyzed in Sec. IIIII For longer 
pulses the distributions change their character. The first 
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(a) 



FIG. 4: Final distribution of the parallel momenta for differ- 
ent initial energy, E and different pulse length (i.e. number of 
cycles in the pulse, n) in the double ionization of the nitrogen 
molecule parallel to the field axis. F(t) = 0.07 a.u., u) = 0.057 
a.u., d = 2.067 a.u. and 9 = 0. (a) E = -0.3 a.u., n = 2; 
(b) E = -0.3 a.u., n = 26; (c) E = -0.6 a.u., n = 2 and (d) 
E = -0.6 a.u., n = 26. 



and fourth quadrants of the panels become strongly pop- 
ulated implying that a number of sequential decays sig- 
nificantly increases. The reason for that is quite obvi- 
ous. After re-scattering when a highly excited two elec- 
tron complex is created, there are two dominant scenarios 
within the first half cycle of the field: a non-sequential 
double escape or a single ionization. During the next 
cycles, unless the molecule is already doubly ionized, we 
are left with a singly ionized molecule which may survive 
to the end of the pulse or the second electron escapes 
and that corresponds to a sequential double ionization. 
The longer pulse duration the more sequential ionizations 
which may easily overcome the number of non-sequential 
events. Thus for sufficiently long pulses, even though 
the re-scattering scenario is involved in the double ion- 
ization process, the momentum distributions will show 
signatures of the sequential electron escape. If the ini- 
tial energy E is much higher than the minimal energy 
of a saddle (Vs ~ —1.2 a.u.) then probability of non- 
sequential ionization is bigger and even if the pulse dura- 
tion is quite long (e.g. n = 26 in panel (b) of Fig. 0} the 
signatures of non-sequential process remain (contrary to 
panel (d) where for E = —0.6 a.u. such signatures are 
not visible). In Fig.[3]one can find the distributions of ion 
parallel momenta that correspond to the data presented 
in Fig. 0] The distributions, as expected, are much nar- 
rower in the case when sequential ionization dominates 
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FIG. 5: Final distribution of ion parallel momenta corre- 
sponding to Fig.gl F(t) = 0.07 a.u., w = 0.057 a.u., d = 2.067 
a.u. and 8 = 0. (a) E = —0.3 a.u., n = 2 - solid line and 
n = 26 - broken line; (b) E = —0.6 a.u., n — 2 - solid line and 
n — 26 - broken line. 



than in the case when non-sequential process does. 

When we change the orientation of the molecule the 
distributions do not change significantly as one can see 
comparing Fig. and Fig. 21 For the field amplitude 
used in the simulations (which corresponds also to the 
experiment [l4|) the saddles for non-sequential process 
are lying quite away from the nuclei and the positions 
and other parameters of the saddles change slightly only 
with a change of 6 — compare Fig. [3J 

Finally in Fig. [7] we show data for the O2 molecule. As 
expected, from the local analysis presented in Sec lIW for 
the same initial energies E and pulse durations as in the 
case of the N2 molecule observed momentum distribu- 
tions are very similar, compare Fig.0] In Eremina et al., 
experiment [l4| the results for O2 differ from the data for 
N2. The latter shows signatures of simultaneous electron 
escape while for N2 it seems that sequential process dom- 
inates. Our analysis indicates that there is practically no 
difference between N2 and O2 if similar initial conditions 
for highly excited complex are chosen. It strongly sug- 
gests that the differences between the observations in the 
two experiments are due to differences in the early stages 
of the excitation process and the nature of the compound 
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r z1 

FIG. 6: Final distribution of the parallel momenta for the 
nitrogen molecule at non-zero angle to the field axis. 
F(t) = 0.07 a.u., u> = 0.057 a.u., d = 2.067 a.u., E = -0.3 
a.u. and n = 2. (a) 9 = ir/4 and (b) 9 = 7r/2 

state before the final decay towards multiple ionization. 

Our analysis also indicates that to increase the ratio of 
non-sequential to sequential ionizations very short pulses 
should be used. It opens a possibility for observation of 
clear signatures of simultaneous electron escape for O2 
molecule, too. 

V. SUMMARY 

To summarize, we have performed a purely classical 
analysis of the final stage of the non-sequential double 
ionization of molecules in the strong laser field. It is 
based on the fact that all trajectories leading to the non- 
sequential escape of electrons have to pass close to the 
saddles in the potential that is formed when the laser 
pulse arrives |a, Ifl L3. l2| • 

We have started with the local analysis of the potential 
within the adiabatic approximation, and we have iden- 
tified and described the saddles. Thereafter, we have 
shown results of numerical simulations. The later allow 
us to draw two conclusions: i) From the point of view of 
classical analysis there is no difference between nitrogen 



and oxygen molecules in a sense that both of them can 
show signatures of simultaneous double escape, ii) Ori- 




FIG. 7: Final distribution of the parallel momenta for the 
oxygen molecule parallel to the field axis for different pulse 
length. F(t) = 0.07 a.u., to = 0.057 a.u., d = 2.28 a.u., 
E = -0.3 a.u. and 9 = 0. (a) n = 2 and (b) n = 26. 



cntation of the molecule with respect to the field axis 
does not influence significantly the final momentum dis- 
tribution for the initial energy range considered in the 
paper. 

The numerical results and their interpretation suggests 
strongly that shorter laser pulses should lead to an in- 
crease of the ratio of non-sequential to sequential ioniza- 
tions. That suggestion could be tested experimentally. 
Higher relative efficiency of non-sequential process would 
be visible in the momentum distributions as a more pro- 
nounced symmetrical escape of the electrons. 
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